Quasiparticle Electronic structure of Copper in the GW approximation 
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We show that the results of photoemission and inverse photoemission experiments on bulk copper 
can be quantitatively described within band-structure theory, with no evidence of effects beyond 
the single-quasiparticle approximation. The well known discrepancies between the experimental 
bandstructure and the Kohn-Sham eigenvalues of Density Functional Theory are almost completely 
corrected by self-energy effects. Exchange-correlation contributions to the self-energy arising from 
3s and 3p core levels are shown to be crucial. 
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Experimental techniques for the determination of the 
electronic bandstructure of solids have made consider- 
able progress in recent years [1]. From the theoretical 
point of view, state-of-the-art calculations within many- 
body perturbation theory allow to obtain band energies 
in a rigorous way, i.e. as the poles of the one-particle 
Green's function G [2]. To obtain G, one needs the elec- 
tron self-energy S, usually evaluated according to the so- 
called GW approximation, derived by Hedin in 1965 [3], 
which is based on an expansion in terms of the dynami- 
cally screened Coulomb interaction. However, due to the 
high complexity and large computational requirements 
of ab-initio calculations of S, the experimental band- 
structures are often compared with the results of (sim- 
pler) calculations performed within Density Functional 
Theory (DFT) [4], in the Local Density Approximation 
(LDA) [5] or in the Generalized Gradient Approximation 
(GGA) [6]. The consequences of this approach (which 
can be set on a firm ground by considering the exchange- 
correlation potential of the LDA or GGA as an approx- 
imation to the self-energy operator) must however be 
considered with great care, particularly when the sys- 
tem under study differs from those -semiconductors and 
insulators- for which the approximation is usually made. 
A recent example is given by Ref. [1], where the mea- 
sured Cu bandstructure is compared with state-of-the- 
art DFT-GGA results, finding discrepancies which are 
both large and dependent on the considered band and 
k-point. As a consequence of this result, the authors of 
Ref. [1] state that it is intriguing to find such pronounced 
"deviations from the band theory" , copper being much 
less correlated than metals with an open 3d shell (e.g., 
Ni). 

In the present work, we demonstrate that copper does 
not, however, deviate from quasiparticle band theory. It 
only deviates from the simple "rigid shift" behavior of the 
self-energy commonly found in the case of semiconduc- 
tors and insulators. For these, the application of the GW 
method to compute self-energy corrections ontop of ab- 
initio DFT results has become a quite well-established 



and standard technique, giving energy levels generally 
in good agreement with experiments, even for compli- 
cated systems like reconstructed surfaces and clusters 
[2,7]. The gaps between empty and filled states generally 
increase by a substantial amount with respect to those 
obtained in the Kohn-Sham (KS) formulation of DFT, 
reaching agreement with experimental results. The KS- 
DFT underestimation of the experimental gap is gener- 
ally found to be weakly dependent on the particular band 
or k-point, although some semiconductor surface gaps 
show stronger QP corrections than bulk ones [9] . This is 
the basis for the introduction of the so-called "scissors 
operator" , often invoked in order to correct the discrep- 
ancies by rigidly shifting upwards the DFT empty bands, 
hence avoiding explicit self-energy calculations. A priori, 
there is no reason why the "scissors" approach should also 
work for metals, where no gap exists between filled and 
empty states. 

Unlike semiconductors, the case of metals has received 
limited attention so far. The band width (i.e., the en- 
ergy range of filled states) of simple metals has often been 
compared with that calculated for the homogeneous elec- 
tron gas (jcllium), finding discrepancies of a few tenths 
of eV [8]. On the other hand, the band structure of tran- 
sition and noble metals cannot be approximated by that 
of jellium: valence electrons of d character play an im- 
portant role both in the electronic structure and in the 
optical properties of these metals. Previous estimates 
of many-body corrections to the 3d Cu bandstructure, 
based on the self-energy of the homogeneous electron 
gas, yielded only partial agreement with the experimen- 
tal data, suggesting that in the case of non-free electron 
metals screening effects (accounted for by the dielectric 
function) should be included in a more realistic way [10]. 
This is precisely the aim of the present work, where a full, 
first-principles GW calculation of the quasiparticle (QP) 
bandstructure of bulk copper is performed. The resulting 
QP energies correct in a non trivial way the DFT-LDA 
results, and are in remarkably good agreement with di- 
rect and inverse photoemission experiments. 
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Among transition metals, full quasiparticle calcula- 
tions have been carried out so far only for Ni [11], while 
we are not aware of similar calculations for noble metals. 
In the case of Ni, GW yields a good description of photoe- 
mission data, except for the 6 eV satellite, which is due 
to strong short-range correlations within the partially- 
filled d-shell. For copper, we expect band-theory to work 
better than for transition metals, since d shells arc com- 
pletely filled. From the computational point of view, 
the very presence of strongly localized d-bands makes 
the ab-initio calculations based on plane waves (PW) 
and norm-conserving pseudopotentials (NCPP) for no- 
ble metals much heavier than those for semiconductors. 

Let us first recall the relations between Kohn-Sham 
eigenvalues, the poles of the one-particle Green's func- 
tion G, and experimental band energies measured in 
photocmission (PE) or inverse photocmission (IPE). In 
DFT, KS eigenvalues cannot be identified with electron 
addition or removal energies, since there is no equiva- 
lent of Koopman's theorem. Instead, the experimental 
bandstructure should be compared with the poles of the 
Green's Function G (r,r';w). The latter are determined 
by an equation of the form: 
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2m ^ r + ^ external ( r ) + VHartree (•*) 



rfr'S(r,r';w) Vnk (r',w) 



ipnk (r, u) + 

E nk (w) V>nk (r, W) , (1) 



containing the non-local, non-hermitian and frequency 
dependent self-energy operator S. The poles of G are 
the QP energies ejjjf , the solutions of ejjjf = E nk (e^f ) • 
Following Ref. [12], we separate the static, bare-exchange 
part S x (r, r') from E c (r, r';w), the energy-dependent 
correlation contribution. Eq. (1) is formally similar to 
Kohn-Sham equations which are solved in the determi- 
nation of the ground-state properties, where the local 
and energy independent exchange-correlation potential 
V xc (r) has been substituted by S. Hence, KS eigenval- 
ues can be considered as a zeroth-order approximation 
to the true QP energies, if the exchange-correlation (xc) 
potential of the DFT is seen as an approximation to the 
true self-energy operator S. This suggests the possibil- 
ity to look for a first-order, perturbative solution of Eq. 
(1) with respect to (S - V xc -lda), where V xc -lda is 
the xc potential of the LDA. We have checked that the 
off-diagonal elements of (£ — V xc -lda) are two orders 
of magnitude smaller than the diagonal ones. To deter- 
minate the QP energies e^f starting from the DFT-LDA 
values we have used an iterative procedure [13]. 

Our bandstructure calculation hence starts with 
a DFT-LDA calculation of the ground-state proper- 
ties, performed using norm-conserving pseudopotentials 
(PPs) and a plane waves basis. The use of soft (Martins- 
Troullier [14]) PPs allows us to work at full convergence 
with a reasonable kinetic energy cutoff (60 Rydbergs if 



the 3s and 3p atomic states are frozen into the core, 
160 Rydbergs when they are explicitly included into the 
valence). Details can be found in Ref. [15], where the 
non-trivial discrepancies between the LDA Kohn-Sham 
eigenvalues and the experimental bandstructure have also 
been evidentiated. In particular, we find substantial dif- 
ferences with respect to the experiment for both the d- 
bands width (3.70 instead of 3.17 eV) and position (more 
than 0.5 eV upshifted in the DFT), in agreement with 
previous calculations [1,16]. 

QP energies are then computed by evaluating S as 
GW [3], so that the dynamically screened potential W{lo) 
(convolution of the inverse dielectric function e _1 (u>) 
with the bare Coulomb potential) is needed. Most GW 
calculations on semiconductor systems use a Plasmon- 
Pole Approximation (PPA) for W(co) [17], based on the 
observation that the Fourier components £q G / (q; u>) of 
the inverse dielectric function are generally peaked func- 
tions of ui, and can be approximated by a single pole. 
Since the evaluation of S c involves an integration over 
the energy, the fine details of the w-dependence are not 
critical, and the PPA turns out to work reasonably well 
for most applications. However, in the case of Cu the use 
of a PPA becomes more critical. The presence of flat el- 
bands 2 eV below the Fermi level implies the presence of 
strong transitions in G , (q; u>) spread over a large en- 
ergy range. For small values of G and G', the behaviour 
of e^Q/ (q; oS) is often very different from that of a single- 
pole function, leading to instabilities when determining 
the Plasmon-Pole parameters. Hence, we have found it 
more convenient to avoid the PPA. Instead, we explic- 
itly compute 6q g , (q; u) over a grid of about 200 fre- 
quencies from zero to ^130 eV, and perform the energy 
integral numerically. Another characteristic of metallic 
systems which leads to additional difficulties in practical 
GW calculations is the discontinuity of the band occu- 
pation crossing the Fermi surface. The problem already 
shows up in the evaluation of the bare exchange matrix 
elements, which read: 

<£f) = (nk|£ x (r,r')|nk) = 



( n k| e i(<i+ G )">ik-q) 2 



q + G 



/ni(k-q)- ( 2 ) 



Here < /„ k < 2 is the occupation of band n at point k, 
and its discontinuity at the Fermi Surface slows down the 
convergency of any numerical evaluation of Eq. (2) on 
a finite q-mesh, no matter how (nk|e^ q+G '' r |nik — q) is 
smooth. To overcome this difficulty, we have first com- 
puted the integral of the term / m (k- q )/ |q + G| 2 over 
small regular volumes centered around each q-point used 
to evaluate Eq. (2), and then used these values in the fi- 
nite q-sum. A similar approach has been used in Ref. [9] . 
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This procedure yields results for (£™ k — V^_ LDA ) com- 
puted with 19 q-points in the IBZ converged within 0.1 
cV, while barely evaluating Eq. (2) over the same set of 
q-points yields an error of about 1.5 eV. 

Normally (e.g., in GW calculations for semiconduc- 
tors) , the calculation of G and W to correct the DFT va- 
lence bandstructure can be performed by including only 
valence states, and fully neglecting the core states which 
have been frozen in the pseudopotential approach. In 
copper, however, when £ is computed neglecting the 3s 
and 3p atomic core states (which in the solid create two 
fiat bands, at about 112 and 70 eV, respectively, below 
the Fermi level), the resulting QP corrections on the d- 
bands are clearly unphysical: GW corrections move the 
highest occupied d-bands above the DFT-LDA Fermi 
level. On the other hand, the situation for s/p states 
(e.g. for the state L 2 ') is much more reasonable, with 
correlation and exchange parts of the self-energy which 
cancel largely each other (as in the case of semiconduc- 
tors), and negative QP corrections of order of the eV. 
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FIG. 1.: Values of (££ k - V" c k ) for bulk copper, plotted as a 
function of the non-interacting energies t n ^. The dotted line rep- 
resents the results obtained without the contributions from the 3s 
and 3p core states. 

The solution of this puzzling situation is provided by 
the role of the above-mentioned 3s and 3p states, which, 
despite being well separated in energy from the 3d ones, 
have a large spatial overlap with the latter. As a con- 
sequence, non-negligible contributions to the self-energy 
are expected from the exchange contributions between 3d 
and 3s/3p states. These contributions can be included in 
our calculation in a straightforward way, by starting with 
a pseudopotential which puts the whole 3 rd atomic shell 
into the valence (see Ref. [15]). Indicating with G^ re 
the non-interacting Green's function containing also the 
3s/3p levels, the self-energy computed as G < Q re W will 
yield the desired dynamical exchange between 3d bands 



and core bands composed of a bare exchange (static) con- 
tribution, E£ ore , plus a correlation part, T, c c ore . Naturally, 
a larger plane- waves cutoff (160 Ry) is needed for conver- 
gence. As a result, we find large bare-exchange contri- 
butions from core levels, leading to very different values 
of (£™ k — Vx^-lda) f° r s /p an d d states, as illustrated 
in Fig. 1. The role of core levels in the calculation of 
the bare exchange contributions (whose importance was 
already addressed for transition metals by Aryasetiawan 
and Gunnarsson [2], but estimated to be of the order of 
1 eV) is hence crucial in the case of copper [18]. More- 
over their effect is unexpected on the basis of DFT-LDA 
calculations, where the band structure does not change 
appreciably even if the 3s/3p orbitals are fully included in 
the valence [15]. The effect of the core states is expected 
to be smaller on the correlation part of S (oj), since the 
screened and bare Coulomb interaction, whose difference 
is involved in ^ ore , are very similar to each other in the 
core region. As a result, Yf° re is small, of the order of 
0.5 eV for d-bands, and close to zero for s/p-bands. The 
inclusion of the correlation part Yf° re does not change 
the widths of d-bands, that remain in good agreement 
with the experimental results, but improves considerably 
their absolute positions [19]. 
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FIG. 2.: Full line: present GW results for the bulk copper band- 
structure, compared with the DFT-LDA results (dashed line), and 
with the experimental data reported in reference [16] (circles). 

In Tab. I we compare our calculations with selected 
experimental data reported in Ref. [16]. In Fig. 2 the 
full theoretical band structure is compared with experi- 
mental data: the agreement is remarkably good and the 
fact that the GW corrections cannot be reproduced by 
any rigid shift of the LDA bands clearly appears. For 
instance, at the L point the shifts change sign for dif- 
ferent valence bands, with QP corrections ranging from 
-0.29 eV to -0.61 eV, to 0.57 eV for the three d-bands. 
As a conseguence, the Li — L 2 > gap decreases after inclu- 
sion of self-energy corrections (see also Tab. I, last line), 
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in contrast with the usual behaviour occurring in semi- 
conductors and insulators. The other gaps at L and X, 
instead, increase. All of them come to good agreement 
with experimental data. Unoccupied (conduction) bands 
at points L, X and K are also obtained in very good agree- 
ment with the experimental data. These findings demon- 
strate that the strong deviations from the single quasi- 
particle band theory in copper, as suggested in Ref. [1], 
do not occur. QP band-theory in the GW approxima- 
tion, instead, remarkably deviates from the DFT-LDA 
+ scissors-operator approach, due to the interplay of the 
different localization and correlation properties of d and 
s/p states. 

In conclusion, we have carried out the first ab-initio 
quasi-particle calculation for a noble metal, copper. In 
order to do that, we had to overcome several prob- 
lems, namely, the large number of plane waves needed 
to describe d orbitals, the calculation of k-space sums 
close to the Fermi surface, and the inclusion of the 
exchange interaction between valence and core levels. 
The single-quasiparticle bandstructure computed within 
Hcdin's GW approximation for the electron self-energy 
turns out to be in very good agreement with experiments. 
The GW method, originally devised to describe the long- 
range charge oscillations [20], is hence shown to yield a 
good description also of copper, a system characterized 
by localized orbitals and short-range correlation effects. 
The corrections to the LDA Kohn-Sham eigenvalues for 
the d-bands are highly non-trivial, since even their sign 
turns out to be k-point and band-dependent, ruling out 
the validity of any "scissors operator" approximation. 
Furthermore, exchange effects between d-electrons and 
3s, 3p core states play a key role and cannot be ne- 
glected in the calculation of QP corrections. Finally, 
a full screening calculation of the self-energy has been 
carried out, avoiding the single-plasmon pole approxima- 
tion. This allows a realistic treatment also of the imagi- 
nary part of the self energy (and therefore the calculation 
of QP lifetimes), which will be studied in a forthcoming 
work. 
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TABLE I. Theoretical band widths and band energies (LDA values + GW corrections, in eV) for copper, at high-symmetry 
points. The striking agreement with the experimental results demonstrates that copper is very well described within a sin- 
gle-quasiparticle band theory, at the GW level. The values in the last column are averages over several experiments, as reported 
in tab. 1 and 13 of ref. [16]. 
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